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A manifestly covariant equation is derived to describe the perturbations 



^ I in a domain wall on a given background spacetime. This generalizes 
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recent work on domain walls in Minkowski space and introduces a 
framework for examining the stability of relativistic bubbles in curved 



Topologically stable field configurations would have been produced during any phase 
transitions that occurred in the early universe. These could be monopoles, cosmic strings 
or domain walls depending on the model. Even a small abundance of these defects can 
have profound cosmological consequences if they were produced during or after inflation. 
In practice, however, calculations are often based on the flimsy assumption that topologi- 
cal defects of high symmetry are the ones which preponderate. Now, perhaps this can be 
justified initially. At least in Minkowski space, the semi-classical approximation to quan- 
tum theory does appear to predict an exponential suppression in the materialization rate 
during a phase transition of classical configurations with low symmetry. [1] However, even 
then, one cannot ignore the fact that there will always be quantum fluctuations attending 
these classical configurations. If one is therefore to place any confidence in cosmological 
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predictions which rely on the assumption that symmetrical topological defects maintain 
their shape in the course of their evolution, such as the collapse of cosmic strings to form 
black holes, one had better be sure no only that the evolution is stable but also that a 
bound can be placed on the perturbation to prevent it from substantially disrupting the 
process. [2] 

In this paper, we examine the evolution of small irregularities in a domain wall prop- 
agating on a curved background spacetime. This provides a generalization of the work 
of Garriga and Vilenkin who examined the stabily of domain walls in Minkowski space 
and equatorial domain walls in de Sitter space. [3,4] Because of the extra technical diffi- 
culty they involve, we defer the examination of lower dimensional defects to a subsequent 
publication. [5] 

We begin in Section II with a derivation of the equations of motion of the unperturbed 
domain wall in the thin wall approximation. We include a discussion of the boundary 
conditions which must be implemented on any physical boundaries. Our approach to 
perturbation theory in Section III will be to expand the action describing domain walls in a 
manifestly covariant way out to second order in the perturbation. We work with the action 
directly rather than with the equations of motion. This not only possesses the advantage 
of preparing the ground for the examination of the quantum theory of fluctuations, it 
also facilitates the identification of appropriate boundary conditions. In Section IV we 
discuss the equations of motion decribing perturbations on various background spacetime 
and domain wall geometries and compare our results with those of Refs. [3] and [4]. We 
work in an arbitrary spacetime dimension, N. 

II. THE EQUATIONS OF MOTION 

Let us consider an oriented domain wall m in the thin wall approximation. This is 
justified so long as the thickness of the wall is much smaller than any other of its dimensions. 
The wall is then described by the timelike hypersurface 
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(2.1) 



p, = 0, • • • , N — 1, a = 0, • • • , N — 2, embedded in spacetime M which we describe by the 
metric g^ u . The metric induced on the world-sheet of the domain wall is then given by 



The action which describes the dynamics of this domain wall is given by the Nambu form 
plus a possible spacetime volume term: 



The first term represents the most simple generally covariant action one can associate with 
the wall, proportional to the area swept out by the world sheet of the wall as it evolves. 
The constant of proportionality a represents the energy density of the wall in its rest frame. 

Before proceeding any further, we need to distinguish between open domain walls 
possessing a boundary, and closed ones which do not. A closed domain wall need not be 
compact. However, if it is not it must be infinite in all directions. A spatial boundary at 
infinity can be ignored. 

An oriented closed wall provides a partition of spacetime into two regions, an interior 
Mi nt and an exterior M ext each supporting its own phase. Neither region need be finite 
in spatial extent. Let p represent the energy density deficit in the region M^. In the 
description of the nucleation of bubbles, M int will be finite. If p is positive (negative), the 
interior consists of true (false) vacuum. We now associate an action with M int proportional 
to the spacetime volume enclosed by the world sheet of the wall. If this volume is infinite 
the associated action will be infinite. However, the change in volume corresponding to a 
variation in the embedding of compact support will always be finite. It will not, therefore, 
affect the equations of motion. 

If, on the other hand, the domain wall possesses a physical spatial boundary it clearly 
cannot provide a partition of spacetime. The phases on either side must coincide and it 
makes no sense to introduce the volume term. 

In the derivation of the equations of motion for the wall, we need to distinguish 
between the physical boundaries dm s of the world-sheet of an open wall and the spacelike 
surface, drat, we introduce to implement the variational principle, and on which the initial 



lab = X^ a X^ b g^ v . 



(2.2) 




(2.3) 
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and final configurations of the domain wall are fixed. One is forced to impose appropriate 
boundary conditions on the spatial boundary. 

The equations of motion of the bubble wall are given by the extrema of S subject to 
variations X — > X^(^) + SX^^) which vanish on dm t : 



5S 



AX» + F^X^X^X^ - pn » = 0, (2.4) 



where A is the scalar Laplacian 

\P 

is the unit normal to the worldsheet and F^g are the spacetime Christoffel symbols 
evaluated on m. We comment on the derivation in the appendix. If p is zero, Eq.(2.4) 
represents a higher dimensional generalization of the geodesic equation describing the 
motion of a point defect. 

Despite the nice analogy, this form of the equations of motion is not very useful in 
practice. This is because all but one linear combination of these equations are identically 
satisfied. To see this, we note that, both on shell and off, 



V6*£ + K^)X%XI = K ab n» , 



where V a is the world-sheet covariant derivative compatible with 7^ and the extrinsic 
curvature tensor K ab is defined by [6] 

K ab = -X%X^, V . (2.5) 

The tangential projections of the Euler-Lagrange derivatives of S therefore vanish identi- 
cally: 

<)S ;X» = 0. (2.6) 



0> ' b 

The geometrical reason for this redundancy is the invariance of the action with respect to 
world sheet diffeomorphisms. In particular, under the infinitesimal world sheet diffeomor- 
phism, 



5S — which implies the (Bianchi) identities, Eq.(2.6). 

It is now clear that the equations describing the world sheet are entirely equivalent to 
the single equation 

oK = p. (2.7) 

If p vanishes, this is just the equation describing an extremal hypersurface. Thus Eq.(13) 
in Ref. [3] which was derived for a domain wall propagating in Minkowski space generalizes, 
with no surprises, to a curved spacetime. That there is only one independent equation 
describing the dymanics of the domain wall is nothing to do with any symmetry, spherical 
or otherwise, that the wall might possess. 

If the wall possesses a boundary, we cannot justify constraining the variation on this 
boundary so that we still have the surface term 

N-2../—c dV =; y ,a 
3m. OA., a 



F-'uV-f^^WXp (2.8) 



to contend with. Here, we have parametrized the boundary £ a = T a (u A ) 1 A = 0, • • • , N — 2 
so that the metric which is induced on dm s is 

fAB ~ du A du Blab - 
The inward pointing normal to d s m in m is — l a . The boundary condition 

X%l a = (2.9) 

is sufficient to ensure that the boundary term vanishes. Physically, this is the requirement 
that no momentum be transferred across the spatial boundary. Modulo the equations of 
motion, this will in turn imply that the boundary of the worldsheet must be a null surface. 



III. THE QUADRATIC ACTION 

At lowest order, the dynamics of any irregularities in the geometry of the wall will 
still be described by the thin wall action Eq.(2.3). One might hope that the non-linearity 
of the equation of motion would serve to damp out any irregularities which might appear 

5 



in the course of the bubble's evolution in much the same way as the non-linearity of 
the underlying field theory is inclined to improve the thin wall approximation in certain 
models. [1] However, the way it turns out (see for example in Refs.[3] and [4]), is that 
sometimes it does but sometimes it does not. 

One way to derive the equation of motion describing the perturbation in the wall is 
simply to consider the linearization of Eq.(2.7) 



with respect to the displacement in the worldsheet, 5X^. This is the method exploited in 
Ref.[3] when the background is Minkowski space. The approach we will follow, will be to 
expand the action out to quadratic order about the classical solution satisfying Eq.(2.7). 
Once this is done, it becomes a simple matter to obtain the corresponding equations of 
motion. In addition, the variational principle provides a guide to the implementation of 
appropriate boundary conditions. 

As we have seen, variations along tangential directions correspond to world sheet 
diffeomorphisms. The only diffeomorphism invariant measure of the perturbation in 
the wall is the scalar 



representing the normal projection of the spacetime displacement vector 5X^. This single 
scalar will now completely characterize the perturbation in the domain wall. 

The simplest way to evaluate the quadratic action is to introduce Gaussian normal 
coordinates for spacetime adapted to the world-sheet hypersurface. Thus, from each space- 
time point P in the neighborhood of m, we drop the geodesic from P which intersects m 
orthogonally at the point P' . P is then uniquely characterized by the coordinates Z, a (P') 
and the proper distance rj along the geodesic. Permitting ourselves an abuse of notation 
which should not lead to confusion, we denote the non-trivial spacetime metric components 



5K = 



(3.1) 



$ = nJX^ 



(3.2) 



by lab(v, £ a )- We note that ^ ab 



(o,e°) = 7«6 



(£ a ) and that 



K ab = -lab 



(3.3) 
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where the prime denotes the proper derivative along the normal and we evaluate it at 
r\ = 0. With respect to these coordinates, $ is simply the component of the variation SX^ 
along the normal, 5X V . The first and second variations in are 

lab (3.4) 
7$ =2K' ab $ 2 - 2<& a <& 6 . 



To second order in $, 



(2) 1 



'-7 = i V-7 



\ ab ^ + l(l {1)2 -^ 1)ab ^ 



(3.5) 



Dropping a divergence, and exploiting the constancy of K implied by the background 
equation of motion, the corresponding second order action can be written 

A^ = --[ d D £^[$A$ + (K' + K 2 )$ 2 ] . (3.6) 

We must however be sure that we can discard the divergence with impunity. This decom- 
poses into surface terms on dm s and dm t . The surface term on dm t causes no problem 
because the appropriate boundary condition in the variational principle is 5<E> = 0. How- 
ever, to ensure the vanishing of the boundary term on dm s we impose the Neumann 
boundary condition 

Z a V a $ = (3.7) 

there. This is the perturbative statement of Eq.(2.9) ensuring that the perturbed surface 
remains null on any finite boundary. 

This form of the second variation is not very useful because it involves K' . However, 
it is simple to express K' in terms of more familiar world-sheet and spacetime scalar s. We 
note that K' appears linearly in the spacetime scalar curvature, N R. The easiest way to 
eliminate K' is therefore to exploit the Ricci identity 

D^D^n* - D v D^n a = N R^pn 13 . 

We contract on [i and a and project onto n v : 

n v D^D v n» - n ■ D(D ■ n) = N R^n u . 
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We now rewrite the first term 



The divergence vanishes because the normal to a hypersurface is a spacetime gradient: 

K' = K ab K ab + N R^v? . 

An alternative way to eliminate K' is to exploit the dynamical Einstein equation for K in 
the initial value formulation of general relativity with the replacement of proper time by 
V-[7] 

We must also expand the enclosed volume out to second order. We find (see appendix, 
Eq.(A2)) 

= - 
2 

We now add Eqs.(3.67) and (3.8) and again exploit the background equation of motion 
to cancel the K 2 term against the volume contribution. 

S( 2 ) = °- f d D t^j [$A$ + { N R^n v + K ab K ab ) $ 2 ] . (3.9) 

2 Jm 

In the elimination of K' we introduced the quadratic in the extrinsic curvature, K ab K ab . 
We can, however, eliminate this term from in favor of curvature scalars by exploiting 
the contracted Gauss-Codazzi equation: 

N R af 3^h a ^ = JV " 1 R + K ab K ab - K 2 , (3.10) 

where the projection, h al3 = g 01 ^ — n a n^. We then get 

S {2) = ^- I d%/ = 7[$A$ + { N R„ v n*n v + N R a fr v h a *tf v -"^R + K 2 )® 2 ] . (3.11) 

2 Jm 

Finally, we use Eq.(2.7) to eliminate K in favor of p and a and add the spacetime curvature 
terms 

N R^n v + N R a ^ u h^h^ = N R^ V , 
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to obtain 



(3.12) 



IV. THE LINEARIZED EQUATIONS 

The equation of motion for small perturbations is now given by 

A$+ ^ N R^h^- N ~ 1 R + $ = 0. (4.1) 

This is a scalar wave equation for $ on the curved background geometry of the world- 
sheet. If this geometry is flat, N R^ va f3 = and we reproduce Eq.(22) of Ref.[3]. The wave 
equation then involves only the intrinsic geometry of the worldsheet. 

Eq.(4.1) is an unconventional Klein-Gordon equation in several respects. 

1. We first note that the perturbation possesses a tachyonic mass whenever p ^ 0. When 
p = this mass is zero. 

2. There is a no n- minimal coupling of the perturbation to the scalar curvature of the 
background world-sheet geometry. This coupling is universal and independent of the 
dimension of the geometry. In particular, there is no privileged dimension in which 
the coupling becomes conformal. 

3. When the ambient spacetime geometry is not flat, the perturbation couples to the 
tangential projection of the Ricci tensor. This is the only explicit dependence of $ 
on the spacetime geometry. It is this feature which distinguishes the perturbation 
theory we are considering from a conventional field theory on the curved spacetime 
described by the metric 7 a &. If, however, the background is a vacuum solution to the 
Einstein equations, this term vanishes. In particular, it will vanish on Schwarzschild 
spacetime. We note also that the perturbation does not couple to the Weyl part of 
the background curvature. 
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Normally, we would interpret a tachyonic mass to signal an instability. However, what 
is more significant is the effective mass given by 

m 2 = N-1 R _N _ ^ 2 (4 2) 

which might be positive depending on the value of the the first two terms. This could 
depend on the topology of the domain wall about which we are perturbing. However, 
even a tachyonic effective does not always signal an instability. This will be the case for 
perturbations about an equatorial bubble in de Sitter space. [4] The expansion of de Sitter 
space introduces a damping term into the Laplacian which annuls the destabilizing effect 
of a tachyonic effective mass. Indeed, sometimes the notion of stability itself is ambiguous. 
An example is provided in Ref. [3] where the stability of true vacuum bubbles in Minkowski 
space is discussed. The Fitzgerald-Lorentz contraction in the perturbation detected by an 
inertial observer in Minkowski space is sufficiently large to render physically divergent 
perturbations of the wall apparently convergent. 

In a vacuum background geometry with a cosmological constant A (this need not be 
de Sitter space), the coupling to the spacetime Ricci tensor simplifies. The background 
Einstein equations then read 

2A 

N R^ = , (4-3) 

so that Eq.(4.1) reduces to the form 

A* + (2 (^) A -»-iR+ (£) 2 ) $ = . (4.4) 

Suppose, in particular, that the background is de Sitter space. We can express N R in terms 
of the Hubble parameter H 

N R = N(N -1)H 2 . (4.5) 

Let us also suppose that p = 0. We describe de Sitter space by Friedman- Robertson- Walker 
(FRW) closed coordinates 

ds 2 = -dt 2 + H~ 2 cosh 2 (H^dttjj^ , (4.6) 
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where dfl 2 N is the line element on a round N sphere. We now consider a spherically 
symmetric domain wall. In general, if the domain wall is spherically symmetric, its world- 
sheet is isometric to an N — 1 dimensional FRW closed cosmology described by the line 
element 

ds 2 = -dr 2 + a(r) 2 dtt 2 N _ 2 , (4.7) 

where r represents the proper time recorded on a clock moving with the wall. The location 
of the wall at any time r is specified by the polar angle x marking the position the N — 2 
sphere of the wall is embedded on the N — 1 sphere. 

dn 2 N _ x = d X 2 + sin 2 xdtt 2 N _ 2 . 

There are two qualitatively different kinds of trajectory. [2] One of these consists of 
trajectories which begin with zero size at the pole x = and grow to a maximum value 
before recollapsing. The other is a bounce which consists of a bubble originating on the 
equator (x = 7r/2) contracting to a minimum value of x before bouncing back to the 
equator. There is a solution with x = 7r/2 representing a domain wall which spans the 
equator. Such solutions can be interpreted as bubbles which materialize from nothing 
through quantum processes [2] and as such are particularly interesting. The world-sheet 
is now an embedded N — 1 dimensional de Sitter space with the same Hubble paramter 

N ~ 1 R= (N -1)(N -2)H 2 . 
The equation describing small perturbations about this solution is then given by 

A$ + (N - l)H 2 <$> = . (4.8) 

This reproduces the expression obtained in Ref.[3]. We stress, however, that the technique 
used in Ref.[3] to derive Eq.(4.2) depended sensitively on the fact that the embedded 
domain wall spanned the equator. 

Even though the effective mass in Eq.(4.2) is tachyonic, this does not appear to be 
significant. 
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In general, the world-sheet of the bubble will not be an embedded N — 1 dimensional 
de Sitter space. If the bubble is collapsing the scalar curvature of its world-sheet will 
diverge. We recall that the scalar curvature corresponding to the line element (4.7) is 
given by 



where the dot refers to a derivative with respect to r. Both a and o remain finite as a — > 
At some point, therefore, the effective mass of the perturbations about the collapsing 
bubble will be rendered real. It is also true, however, that as we approach the final stages 
of collapse, the thin wall approximation will break down. 



We have provided a framework for the examination of perturbations of domain walls 
on a given spacetime background. 

This analysis can be extended in at least two different directions. 

The first is the treatment of perturbations on lower dimensional topological defects. [5] 
When the co-dimension of the world-sheet is r, there will be r scalar fields describing the 
perturbation, one for the projection of onto each of the r normal vectors M . What 
is more, the equations we obtain will generally be coupled in a non-trivial way. To derive 
the equations of motion, we need to develop a different line of attack. The reason for this 
is that for co-dimension r > 1, we can no longer exploit a Gaussian system of coordinates. 
The description of the extrinsic geometry will, in general, be more complicated. Now there 
will be r extrinsic curvature tensors, one for each n^ M : 



which vanish on a hypersurface now appear with a vengeance. [6] 

The weak point in our treatment of perturbations is that the domain wall has not 
been treated as a source for gravity. This is a serious limitation. If we are to place any 




V. CONCLUSIONS 



In addition, so-called torsion terms of the form 
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confidence in perturbation theory, we need to accommodate the back-reaction. When this 
is done, the simple extremal form Eq.(2.7) gets replaced by the Lanczos equations, 

AK ab = 8nGa<y ab , (5.1) 

where Aif a (, is the discontinuity suffered by K ab across the domain wall. [8] These equations 
are very different from Eq.(2.7). If they do possess Eq.(2.7) as their limit when the coupling 
to gravity is turned off, this is not obvious. What is more, whereas the solution of Eq.(5.1) is 
relatively straightforward when the domain wall is spherically symmetric[9], the treatment 
of perturbations about such a wall is far from trivial. For now, the displacement in the 
wall 5X^ will couple to perturbations in the spacetime metric with the generation of 
gravitational waves. This is currently being examined. 
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APPENDIX 



The extrema of S with respect to variations which vanish on the boundary satisfy the 
Euler Lagrange equations 



Now 



so that 



and 



d <9^/~7 d\f~l 



di a dX% dX» 



ab 



dial 



~dXj' 



= 77 



^7° , 



(Al) 
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We also have 



dXv 2 

The first derivatives of the spacetime metric appear in the combination T^. The term 
in square brackets Eq.(Al) reproduces the corresponding term appearing in Eq.(2.4). To 
complete the derivation of Eq.(2.4), we note that under the variation — > X^ + 5X^, 
the volume transforms by 



5V 



or alternatively, in the notation of section III: 

5V = J d N -^y/^y^ . (A2) 
In this form, it is clear that the second variation is given by Eq.(3.7). 
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